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We demonstrate that the MSSM flat directions can naturally account for the seed magnetic fields 
in the early Universe. The non-zero vacuum expectation value of an MSSM fiat direction condensate 
provides masses to the gauge fields and thereby breaks conformal invariance. During inflation the 
condensate receives spatial perturbations and SU(2) x U(1)y gauge currents are generated together 
with (hyper)magnetic fields. When these long wavelength vector perturbations reenter our horizon 
they give rise to (7(l) cm magnetic fields with an amplitude of 10 -30 Gauss, as required by the 
dynamo mechanism. 
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A large scale magnetic field with a coherence scale of 
few kpc and an amplitude of order \iG is observed today 
in galaxies and galactic clusters . The origin of the cos- 
mic magnetic field may well be found in the small primor- 
dial magnetic fluctuations in the early Universe, which 
are amplified by the galactic dynamo mechanism P, | 
While various astro-physical processes, such as battery 
or vorticity effects |^ face an uphill task in sustaining the 
large scale coherence length j^J, there is no shortage of 
microphysical processes in the early Universe that could 
generate a magnetic field . The problem is usually the 
shortness of the magnetic coherence length, which how- 
ever could be enhanced by the inverse cascade of magne- 
tohydrodynamic turbulence Q. 

It has been noted already many years ago that the 
coherence length problem might be solved by inflation 
provided conformal invariance gets broken during infla- 
tion. Interesting possibilities have been explored where 
the conformal invariance is broken either via triangle- 
anomaly of the gauge field stress energy tensor ||, or 
via the gravitational or anomalous coupling of the pho- 
ton to the axion 9]. There are also proposals for excit- 
ing the massive Z-boson of the Standard Model during 
inflation and preheating, which could produce [/(l) em 
magnetic field after the electroweak phase transition via 
hypercharge mixing |To| . 

In this paper we will discuss a simple mechanism 
for generating a seed magnetic field during inflation 
in the context of the Minimally Supersymmetric Stan- 
dard Model (MSSM). As is well known, MSSM has 
flat directions, made up of gauge invariant combina- 
tions of squarks and sleptons, which may acquire non- 
vanishing vacuum expectation values (vevs) during infla- 
tion, thereby forming homogeneous zero-mode conden- 
sates. The condensates may play a significant role in 
many cosmological phenomena [llj, such as generating 
baryons and dark matter particles by first fragmenting 
into Q-balls which then decay [J^j ^ has also been sug- 
gested that the origin of all matter and density perturba- 
tions could, in principle, be due to MSSM flat directions 
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By definition, flat directions are the minimum energy 
configurations even when supersymmetry is broken. Be- 
cause of the condensate the SM gauge fields obtain non- 
vanishing mass terms, thereby breaking the conformal in- 
variance. This is not by itself sufficient for generation of 
magnetic fields, since for the homogenous condensate the 
gauge field configurations are pure gauge However, 
during inflation the flat directions obtain quantum fluc- 
tuations from the inflaton induced false vacuum, which 
perturbs the classical condensate and also induces fluc- 
tuations in the gauge degrees of freedom. These may no 
longer be gauged away. The situation may be compared 
with plasma physics, where an electromagnetic field is 
generated in an electrically neutral plasma by virtue of 
the motion of the electrons. Here the condensate pertur- 
bations represent also perturbations in the charge density 
and by virtue of their motions, generate gauge currents. 

While inflation lasts, the gauge field perturbations are 
stretched outside the horizon. When they eventually re- 
enter they provide us with a seed (hyper)magnetic field, 
which after electroweak phase transition is projected onto 
the regular £7(l) em magnetic field. 

Let us consider one particular flat direction with com- 
ponents carrying SU(2) x U(l)y charges such as LLe. 
The field content of this direction is given by 
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where i ^ j , is a complex scalar field and square roots 
are to obtain canonical normalization. The background 
gauge fields are assumed to vanish. With the field con- 
figuration Eq. |Q it is possible to show that the gauge 
invariant matter currents 
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vanish for vanishing gauge fields in a Robertson- Walker 
background g^ = a 2 (r])rj^, where a{rf) is the scale factor 
in comoving coordinates and r\ the conformal time. The 
covariant derivative is given by D p — d p — igA A T A — 
ig'YA^, where A/j,(A A ) is the Abelian (non-Abelian) 
gauge field corresponding to U{l)y (SU(2)). The sum- 
mation in Eq. @ is over the scalar fields belonging to 
the representation of the corresponding gauge group. 
The corresponding equations of motion are given by 



V ^d p F^=a( V ) 2 jY , 
dV 



0. 



(3) 



Here TL is the comoving Hubble expansion rate, 
and the adjoint covariant derivative is determined by 
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F A = d n F, 



P ^ u ^gf ABC A B p FC u , where f ABC is the 
structure constant of the SU(2) gauge group, and F A V is 
the gauge field strength. 

Perturbing Eq. to first order we get 1 



rf"d p (d p SA A - d u 5A A ) = a 2 5J i 
r,P v d p {d p 5A u - d v 6A M ) = a 2 8Jl 
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where the background equation A A = A^ = has been 
used. In calculating the perturbations of the matter cur- 
rents the perturbations of the flat direction, 6(j>, do not 
appear since the current vanishes for all values of <p for 
vanishing gauge fields, so only the terms proportional to 
gauge field perturbations contribute: 



SJ A = ig [L\T A {-igT B 5A B - ig'Y^A^U 

-L\(igT B 5A B + ig'Y L 8A^)T A L t + (L, -> L 

<5jJ = ig'[2Y L Ll(~igT B 6A B ~ig'Y L 8A p )U 

+ L J )-2ig'Y 2 8A tl \e k \ 2 ] (5) 

Inserting the background field configuration of the flat 
direction Eq. into Eq. (JSJ), we finally obtain 

1 



SJ? = ~9'WSA A 
5Jl = g' 2 \4\ 2 6A p 
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Hence Eq. (jHJ) produces effective mass terms for Eq. Q. 
The equations of motion are similar to massive scalar 
field equations of motion and thus in the Coulomb gauge 



SA A = 5A = d i 5A A = d i 8A i = 0, Eq. Q may be 
written as 

(d^-cPd^SAf + ^faivfSAf = 

{dl-Vd^Ai + g'^a^YSAi = 0, (7) 

where i — 1,2, 3. Eq. (JJJ) reminds us of a scalar field 
equation with mass m — \g(f>\/y3 and m! = \g'4>\- The 
solutions can be written as plane wave expansions 



SAf = 



J2 (Mv)e?Wx + w k ( V ye AX *al kx ) e ikx , (8) 



where a corresponding expansion holds for 5Ai, too. Here 
e AX are the polarization vectors which obey k l ei — in 
the Coulomb gauge. The mode functions, Wk{f]), can be 
given in terms of Hankel functions as 
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and similar expression holds for the mode function of SAi 
with v — ► v 1 . In the above expression k is the wavenum- 
ber and 
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where Hj is the Hubble expansion rate during inflation. 
The mode functions 10 reduce to the vacuum mode func- 
tions in the limit rj — > — oo, 
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We are mainly interested in the super Hubble horizon 
fluctuations of SAi,SA A . When k\rj\ <C 1, we obtain, 
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and a similar expression for |<L4j| but with v — » v' . 
The amplitudes of the long wavelength fluctuations are 
frozen and start to re-enter once inflation comes to an 
end. For the sake of simplicity we assume here that 
post-inflationary era does not trigger super-Hubble fluc- 
tuations in the gauge fields which would alter our result 
in Eq. (|12[1. However, note that the vacuum structure is 
different during inflation and the radiation era 2 . There- 
fore at late times it is important to match the above mode 
functions (valid during inflation) to a linear combination 



1 We have neglected the perturbations of metric and non-flat de- 
grees of freedom. In the gauge SA^ = and 9*<5A^ = the 
result for the gauge field equation turns out to be the same. 



2 Note that the flat direction condensate amplitude slides both 
during inflation and post-inflationary era. Therefore the masses 
m, ml vary with time. The MSSM flat direction decays late but 
well before the electroweak scale 1 1 ~J1 . and by decaying they re- 
store the conformal invariance. 



2 



of traveling waves in the radiation dominated epoch via 
Bugoliubov coefficients, a k , 0k, 

w k {r)) = a k w ( ™ c) (Ti) + f31 k w^ ac} (r))* , 
d n w k (v) = a k d v w k vac) (77) + f3*_ k d n w k vac) (77) , (13) 

For simplicity let us assume that the transition from infla- 
tion to radiation happens instantly, and we also assume 
that MSSM flat direction vev becomes zero right after 
the end of inflation so that the condensate disappears 
instantaneously. Although the MSSM flat direction vev 
always goes to zero, it need not vanish right away after 
the end of inflation. However, it is possible to find an ex- 
ample within D-term inflationary models where the flat 
direction vev can become zero right after the end of in- 
flation; this happens due to a large Hubble induced mass 
correction which causes the condensate field to roll down 
to the origin within one Hubble time (lif. 

In this letter we proceed with the assumption of a sud- 
den transition from inflation to radiation on super hori- 
zon scales, which leads to 
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and a similar expression for \(3 k \ 2 with v' — > v. Therefore 
during the radiation era the amplitude of the gauge fields 
is enhanced by a factor (3 k compared to the Minkowski 
vacuum. Note that during radiation era, w k vac ^ (77) oc 
\J\j2ke~ lkTI . Therefore the final the gauge field spec- 
trum ex fc _1_1/ . For v w 1/2, the spectrum is indeed 
flat, similar to the case of massless scalar field. Similar 
conclusions were also reached in Ref. 

During radiation era, at scales where diffusion can be 
neglected, the hypercharge field is frozen in the sense that 
the lines of force move together with the fluid, Ai, Af oc 
a(i)- 2 . 

In our case, the hypercharge field is a linear combi- 
nation of the fields, ~ B^ + Z^. We are mainly 
concerned with the U(l) em photon field, B^, below the 
electroweak scale, given by /3 M = sin^A 3 + cos By/A^, 
where singly ~ 0.23 is the Weinberg's angle. Since the 
non-Abelian field A^ obtains a screening mass in a radi- 
ation dominated Universe [l^ . its contribution vanishes 
and the unscreened photon field is Bi w cos 9\yYi and the 
real magnetic field is Bi = Cij k djB k . 

Let us now estimate the amplitude of the magnetic 
field at the length scale I where k = 2ir/l <C Hj. We 
obtain 



Bi « cos 6^ CV) 
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where 
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Within our setup LLe obtains a vanishing vev right after 
inflation so that the amplitude of the magnetic field is 
frozen during the radiation epoch, evolving as B ~ a~ 2 . 
Since the temperature scales as the inverse of the of the 
scale factor, the magnetic field at the time of galaxy for- 
mation is given by Bi^ — Bi(T g f /Tr) 2 , where T g f is the 
temperature at the galaxy formation and Tr is the reheat 
temperature. The length scale also evolves, / ~ a, so the 
length scale at the galaxy formation is related to the one 
at the end of inflation by l g f = l{Tn/T g f). Putting every- 
thing together we obtain the amplitude of the magnetic 
field at the time of galaxy formation as 

B i,gf ~ cosO w C{v)\--v \Hj< l gf 



Til 
Tr 



2 

1/2+1/' 



(17) 



We estimate the reheat temperature by assuming that 
the inflaton potential energy density is directly converted 
to a thermal bath, V end w 3M 2 H] « (7r 2 /30)g >t T^, where 
,g* = 106.75, is the number of relativistic degrees of free- 
dom. Thus we obtain an order of magnitude estimate for 
the magnetic field which reads 



B, 



gf 



7- 10" 30 G (l + zgf)" 1 
T \ fl kpc 
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2.73 K 
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when (1/2 - v') ra (g'^/Hj) 2 < 1/2 and T gf = (1 + 
z g f)T , where z g t is the redshift of the galaxy formation 
time and a factor of 5 • 10 4 (1 + z g f)~ 2 due to the 
increase of flux at the collapse of protogalaxy is included. 

Note that the the strength of the magnetic field de- 
pends on the vev of the hypermangetic field during in- 
flation. Obviously when the vev is zero, the conformal 
invariancc is restored and there is no magnetic field. For 
a field almost massless during inflation, we would natu- 
rally expect that the MSSM flat direction vev is locked 
to the expansion of the Universe such that \<j>\ ~ Hj. 
This yields the maximum amplitude of the primordial 
magnetic field to be of order 10~ 30 G with a coherence 
length I ~ 1 kpc as required in [lflj . 

A word of caution, though: the above estimation of 
the strength does not take into account of the sub- 
horizon plasma physics due to charge separation and in- 
homogeneous current distribution, see for instance |20| . 
Note that MSSM flat directions do not generate cur- 
rents at any time but nevertheless the current density 
obtains non- vanishing fluctuations, which requires a de- 
tailed analysis of the sub-horizon effects which goes be- 
yond the scope of the present paper. On the other hand 
the condensate is assumed to decay instantly at the end 
of inflation, so that conformal invariance of the gauge 
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field is restored. In reality one would expect the conden- 
sate to be stuck at its inflationary value until the Hubble 
parameter is of the order of its mass ~ ITeV. Thus the 
magnetic field amplitude decay rate is slower than in the 
conformally invariant case. Also in reality the condensate 
should be slow-rolling during and after inflation unless it 
is at the minimum of the potential, but we expect that 
this produces only minor quantitative changes to our re- 
sults. 

In conclusion, we have provided a natural mechanism 
within supersymmetry for generating seed magnetic field 
based on generating non-vanishing masses to the SM 
gauge fields from the vev of the MSSM flat directions. 
The simplicity of our scenario lies in that we do not break 
conformal invariance by an ad-hoc assumption, nor do we 
invoke any exotic fields or coupling other than the ones 
present in the minimal extension of the SM. The excited 
gauge fields, mainly the hypercharge field arising from 
LLe flat direction, are stretched during inflation outside 
the horizon in a similar way as the scalar fluctuations 
responsible for generating the perturbations in the cos- 
mic microwave background radiation. The subtle point 
here is that we generate non-vanishing current density 
required for breaking the conformal invariance through 
the fluctuating current. 
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